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Abstract 

We extract the neutron charge form factor from the charge form factor of deuteron 
obtained from T2q{Q 2 ) data at < Q 2 < 1.717 (GeV 2 ). The extraction is based on the 
relativistic impulse approximation in the instant form of the relativistic Hamiltonian 
dynamics. Our results (12 new points) are compatible with existing values of the 
neutron charge form factor of other authors. We propose a fit for the whole set (35 
points) taking into account the data for the slope of the form factor at Q 2 = 0. 
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The behavior of the neutron charge form factor G^Q 2 ) , (Q 2 = — q 2 , q - the momentum 
transfer) is of great importance for the understanding of the electromagnetic structure of 
nucleons and nuclei. However, G%{Q 2 ) is still known rather poorly. 

As there are no free neutron targets, G%{Q 2 ) has to be extracted from the data for 
composite nuclei, for example deuteron or 3 He |, g, |, |, g, |, % |, % [TO], 0, [JJ, . The 
direct measurement of great precision (~ 1.5%) is possible only for the slope dG%{Q 2 )/dQ 2 
at Q 2 = 0, as determined by thermal neutron scattering JT1 . 
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While obtaining the information about the neutron from the scattering data on compos- 
ite systems one encounters two kinds of difficulties. First, the results depend crucially on 
the model for NN interaction |T3|, ITR [TB| . Second, there exists a dependence on the rel- 
ativistic effects, exchange currents, nucleon isobar states, final state interaction in inelastic 
channels etc. ML 0|. The use of polarized beams and polarized targets in recent experiments 



diminishes uncertainties due to those effects |3|, B, |^, ^, H, |], |TT|, [T3| ]. 

In the present paper the neutron charge form factor is extracted from the experimental 
data on the deuteron charge form factor obtained through polarization experiments on elastic 
ed scattering [18], B5|. In the JLab experiments Jl9| the deuteron charge form factor is 
obtained up to Q 2 = 1.717 (GeV 2 ). In this range of momentum transfer the theoretical 
description of the polarization tensor T2o(Q 2 ) depends essentially on the choice of the form 
of NN interaction and relativistic approach is required. 

Our calculations are based on the method of relativistic Hamiltonian dynamics (RHD) 
which is widely used in present time. One can find the description of RHD method in 
the reviews [[2t]] and especially the case of the deuteron in the reviews |2T 
use our own variant 



23], |24] of the instant form of RHD. This variant permits to take 
correctly into account the relativistic effects in the elastic ed scattering in the relativistic 
impulse approximation P5| . The main point of our approach is the construction of the 
electromagnetic-current operator for the system of interacting particles. In our approach 
this operator is Lorentz covariant and satisfies the conservation law. 

Let us note that, as far as we know, it is for the first time that the neutron charge form 
factor is determined from an analysis of the deuteron charge form factor. 

In our approach in the relativistic impulse approximation the following equation for the 
deuteron charge form factor takes place (see [2"5] for details): 



We 



G C m(Q 2 ) G p m (Q 2 ) + G n M (Q 2 ) 



V, 



G C (Q 2 ) = Gcc(Q 2 ) [G P E (Q 2 ) + G n E (Q 2 

Here G p E n M are charge and magnetic form factors of proton and neutron. The fact that 
nucleons magnetic form factors enter the Eq. ([I]) is due to the relativistic effect. 
The functions Gcc , Gqm m (0) are given by: 

Gcc(Q 2 
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here <p l (s) is the wave function in the sense of RHD (see 
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(s) = tfsui(k)k , k = -Vs-AM 2 
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M is the nucleoli mass, / = 0,2 - the nucleoli angular momentum in the deuteron, ui(k) - 
the wave function for the model NN interaction. The functions g l Q C , g l Q M are given by the 
following equations (^-([TT]) (note that the same equations were obtained independently in 
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L 2 = L 2 (s, Q 2 , s') = P21 (Pi 2 - 6 - P21 (P22 - 6 P 
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Here &% and u>2 are the Wigner spin rotation parameters: 

£(s,QV) 



U\ = arctan ■ 



co>2 = arctan ■ 



» r {Vs + v^ 7 ) 2 + Q 2 \ + + v^ 7 ) 

a(s,s%(s,Q 2 ,s') 



M(s + s' + Q 2 )a(s, s') + Vss'(AM 2 + Q 2 ) 

and a(s, s') = 2M + y/s + v 7 ^ 7 - 

P2« = p2i(z) , P'u = p2i(z') , i = 0, 1, 2 - the Legendre functions: 

P»W = i(3^ 2 -l) , P 21 (z) = 3zv / I^ 2 , P 22 (z) = 3(l 



(7) 
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z = z(s, Q , s') 
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z' = z'(s,Q 2 ,s') = -z(s',Q 2 ,s). (9) 
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tf(s, Q , s') = 6(s' — s\) — 6(s' — s 2 ), is the step function. 



Q 2 (Q 2 + 4M 2 )s(s - 4M 2 ) . (10) 



The functions Si^s, Q 2 ) give the kinematically available region in the plane (s, s') (see [27]). 



g l < l (s,Q\s')=gy t (s',Q\s 
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Using Eq.(|l]) one can write the neutron charge form factor in the form: 
G C (Q 2 ) G CM (Q 2 ) 



Ge(Q ) 



Gcc{Q 2 ) Gcc{Q 2 



G P M (Q ) + G n M (Q 2 ) — G P E (Q ) 



(12) 



We calculate the nucleon charge form factor in the points Q 2 where the deuteron charge 
form factor GciQ 2 ) is measured. In these points the nucleon form factors G P E (Q 2 ) , G P M {Q 2 ) , 
Gm(Q 2 ) are obtained through the fits of their experimental values. The functions GcciQ 2 ) , 
Gcm{Q 2 ) can be calculated using the equations @, (§) and some deuteron wave functions. 

Let us discuss now the problem of choosing the deuteron wave functions to use for the 
calculation of GcciQ 2 ) > Gcm(Q 2 ) QT^D- We investigated the behavior of T 20 (Q 2 ) (our 



determination of T2o(Q 2 ) is the same as in [0) using different wave functions]]: Paris wave 
functions [p8| |, the versions I, II and 93 of the Nijmegen model [29|, charge-dependent version 
of Bonn potential [BO], and the wave functions obtained by the relativistic dispersion variant 
of the inverse scattering method [31]. As T 2 o(Q 2 ) for polarized ed scattering depends weakly 
on the form of nucleon form factors one can use the experimental data for T 2 o(Q 2 ) to choose 
the most adequate deuteron wave functions. Fig.l presents the results of our calculation of 
T2o(Q 2 ) with the use of the wave functions 0, |29], |30], 

|i~5| as well as the experimental points from the papers |17l [18], [19], [32], [33], 0, [35] 



and nucleon form factors from 



One can see that the best description of T 2 q(Q 2 ) is given by the wave functions pTT 



Our estimations show that other wave functions (e.g. used in |36| , 
description of T2o(Q 2 ) than the wave functions [31 



(]]) also give poorer 



Let us emphasize that the wave functions [[H|] were obtained more than 20 years ago and 
so no possible fitting reasons for T 2 o{Q 2 ) could influence the choose. These wave functions 
used in the relativistic calculation of the function A(Q 2 ) give the correct behavior up to 
Q 2 ~ 3 (GeV 2 ). 



It seems to us that the validity of the wave functions [[JT|] is due to the fact that they 
are "almost model independent": no form of iVW interaction Hamiltonian is used. The 



1 The details will be published elsewhere. 
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wave functions |3T| were obtained in the frame of the potentialless approach to the inverse 
scattering problem (see for the details ||38|| ). They are given by the dispersion type integral 
directly in terms of the experimental scattering phases and the mixing parameter for NN 
scattering in the 3 Si — 3 D\ channel. 



one 



In the Eq.([L2|) we use for the nucleon form factors G P E {Q 2 ) , G P M (Q 2 
(with the best x 2 ) of the fits of the recent paper f3pj] - DRN-GK(3). 

The results of our calculations of the neutron charge form factor in the points where the 
deuteron charge form factor is measured are given in the Table 1 (see also Fig. 2). 

The accuracy of our calculations are determined by the accuracy of measurements of 
charge deuteron form factor |17|, [18], |19| and nucleon form factors which are the folowing at 
Q 2 < 1.717 (GeV 2 ): for G P E (Q 2 ) 1-10% § H H @, for G P M (Q 2 ) 1-3% 
for G n M {Q 2 ) 1-10%. |, ^ [§, 0, g 



\ 42, 43, 44 



We obtain the first three points at low momentum transfer from the data for the deuteron 
charge form factor given in the paper |1| . In this range of momentum transfer the behavior 
of the deuteron charge form factor and so G E (Q 2 ) do not depend on the choose of the wave 



functions [28, 2£,m, 31 



The first point at Q 2 ~ 0.16 (GeV 2 ) is almost the same as in ||, however, our errors 
are much smaller. The second and the third points are compatible (within the experimental 
errors) with the points of |7], [|, |Il|] . Our point # 7 is in fact the same as in [|IJ but our error 
is larger. 

Our values of G E in other points (at Q 2 > 1 (GeV 2 ) are strictly positive. This result 
differs from e.g. the results of the paper Q consistent with G% = 0. Let us note that our 
errors at Q 2 > 1 (GeV 2 ) are sufficiently small, smaller than, e.g. in |], |J. 

Our values # 4-8 are extracted from the values of charge deuteron form factor of the two 
different works |17], 19]. The results of these works are in rather poor agreement with each 
other in the region of the first dip. So the values of # 4-8 of G E are not well determined in 
the present work. One needs additional experiments in this region. 

It is now interesting to fit all the existing values of neutron charge form factor ([p], [|. || 
H |], [7], H, §, pi], [ll]] and Table 1). We use for the fitting the following function (see Jl5| and 
the review [2~T]) with two parameters a and b: 



G n E {Q 2 



a t 



Air, 



l + br 



G D {Q 2 ) , G D {Q 2 



1 + 



0.71 



Q' 2 



AM 2 



(13) 



The neutron magnetic moment //„ = -1.91304270(5) []|J. Q 2 in G D {Q 2 ) is given in (GeV 2 ). 
We obtain the parameter a from the slope of the neutron charge form factor at Q 2 = 
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CI 1391: 

Ann 



dQ 2 



= 0.0199 ± 0.0003 fm 2 . (14) 

=o 

The fitting of the slope ( |14D gives a=0.942 with the accuracy ~ 1.5%. 

This value of a gives the slope of G^(Q 2 ) at Q 2 = which is measured directly in the 
experiment. 

The parameter b is fitted using the x 2 criterion. If we use all the 35 points we obtain b 
= 4.61 with x 2 = 69.0. Note that the fit DRN-GK(3) g§ of 23 points has x 2 = 63.9. 

If we exclude the points # 4-8 then the 30-point fitting gives b = 4.62 with x 2 — 61.5. 
As the errors of these points are large this fitting differs from the previous one slightly. 

Let us note that our fitting for 23 points of the papers P], |3|, [|, [^, [], [I], |, [| |T(], [□]] (not 
taking into account our points) gives b = 4.69 with x 2 — 57.7. The two curves lie near one 
another (see Fig. 2) so our points are consistent with the known points of other authors. 

The results of fitting, the experimental points |1], |, ||, |6], 0, || ^, |TD|, 0, as well as our 
new points are shown on the Fig. 2. The points # 5 and # 6 are out of the figure. 

To summarize, 

1) We extract new points for the neutron charge form factor from the experimental data 
for the deuteron charge form factor. The obtained values are consistent with the known 
values of other authors. 

2) We perform the fitting for 35 values of the neutron charge form factor including our 
points. The fit has the form fllTf ) with a = 0.942, b = 4.61. 

This work was supported in part by the Program "Russian Universities - Basic Re- 
searches" (grant # 02.01.28). 
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Table 1. The values of G^(Q 2 ) / G d(Q 2 ) obtained in the present paper. The values of the 
deuteron charge form factor used for the extraction of G^(Q 2 ) are also given. 



# of points 


Q 2 (GeV 2 ) 


Gc(Q 2 ) 


Ref. 


G%{Q 2 )/G D {Q 2 ) 


1 


0.160 


0.163±0.017 




18 




0.076±0.116 


2 


0.215 


0.100±0.012 




18 




0.052±0.129 


3 


0.303 


0.035±0.020 




18 




-0.234±0.401 


4 


0.556 


(0.127tH56 7 )"10- 1 




17| 


1.23±0.92 


5 


0.651 


(-0.117±0.162)-10- 2 




19 




-2.61±1.65 


6 


0.693 


(0.166™)-l(r 2 




17 




-4.52±5.10 


7 


0.775 


(-0.253±0.063)-10- 2 


f 


19 




0.677±0.361 


8 


0.831 


(-0.147i^)-10- 2 


I 


17 




-0.140±0.432 


9 


1.009 


(-0.396±0.028)-10- 2 




19 




0.389±0.107 


10 


1.165 


(-0.348±0.031)-10- 2 




19 




0.259±0.131 


11 


1.473 


(-0.310t°;^)-10- 2 




19 




0.405±0.263 


12 


1.717 


(-0.194^1).10- 2 








0.174±0.294 
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Figures capture. 



Fig.l. Data and the results of calculation of the deuteron polarization tensor T2o{Q 2 ) for 
the elastic ed - scattering with the use of the nucleon form factors from the paper |D| and 
different wave functions. The experimental points are: open circles - ||32[| , open squares - 
351, open triangles - p4], filled circles - [O], filled squares - [fLSfl, filled diamonds - |19|, filled 



triangles - |p3 |. Curves: solid - Nijmegen-II [ 29] , dashed - [[U|], dotted - p8 |, dash-dotted 
- Nijmegen-I PSI, dash-dotted-dotted - 



Fig. 2. The experimental values and the results of fitting for the neutron charge form 
factor. The experimental points: bold cross - [Q, open bold diamonds - open up 

triangles - |§, open circles - [|J, open down triangles - ||, open stars - [fLOl , filled circles - 
0, filled diamonds - ||, filled up triangles - [[J, filled stars - filled squares - the present 
work. The points # 5 and # 6 are out of the figure. The curves : solid - the result of fitting 
of 35 experimental points (including our points of the Table 1) using the equation (a = 
0.942, b = 4.61 with x 2 — 69.0), dashed - the result of fitting of 23 points of other authors 
(a = 0.942, b = 4.69 with X 2 = 57.7). 
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FIG.l. 
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FIG.2. 
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